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Abstract 

Using 15.1 pb _1 of data collected by CMD-2 in the c/>-meson energy 
range, the branching ratio of the conversion decay 4> — ► 7r°e + e _ has 
been measured for the first time: 

B(<p -» TT°e + e~) = (1.22 ± 0.34 ± 0.21) • 10~ 5 . 
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1 Introduction 



Conversion decays of a vector meson V into a pseudoscalar meson P and 
a lepton pair l + l~ (V — > Pl + l~, where I — e,/_t) are closely related to the 
corresponding radiative decays of V into P and a photon (V — > P"f). 

In conversion decays a squared invariant mass of a lepton pair 
Mf nv {l + l~) = q 2 or a squared mass of the virtual radiated photon is not 
equal to zero as for usual radiative decays. By studying M inv (l + l~) spectra 
one can determine a so called transition form factor Fp(qf, g|) of pseudoscalar 
mesons P as a function of qf. 

Various theoretical models, from standard vector meson dominance 
(VMD) to calculations on the lattice, predict the transition form factor and 
the resulting branching ratio of the decay [ III El Bl ffl . 

Precise values of the branching ratios of conversion decays are also im- 
portant while studying the yield of direct lepton pairs in heavy ion collisions. 
As it was noted long ago, an observation of the anomalously large yield could 
indicate the existence of the quark-gluon plasma [ ||. In general, both the 
yield of dileptons and their mass spectra are of interest while studying the 
change of the meson properties in medium [||. Recent experiments studying 
dileptons in heavy ion collisions at CERN reported on an excess of the num- 
ber of e + e~ [ [7j and fi + {i~ pairs [ || over the expectations from usual hadron 
decays (it — > e + e~7, 77 — > e + e~7, rf — > e + e~7, u — > ir°e + e~, p/uj —> e + e~, 
<p — > e + e~ as well as similar decays into muon pairs). Its explanation requires 
the accurate calculation of the production rate of pseudoscalar and vector 
mesons as well as good knowledge of both the invariant mass spectra and the 
branching ratios of their decays into lepton pairs. 

The experimental information on such decays is rather scarce [ 0. While 
for the uj meson both possible decays into 7r° were observed, uj 7r°/n + /i _ 
[ and uj — ► ir°e + e~ [ O, only few events of the decay <fi — ► r/e + e _ were 



detected [ 12 



A large data sample of the mesons collected by two detectors at 
VEPP-2M allowed reliable detection of the decay mode — > r]e + e~ by both 



CMD-2 [ p[ and SND [ HJ groups. CMD-2 has also reported on the first 



observation of the — > 7r°e + e - decay [ [13 

This work is devoted to the determination of the branching ratio for the 
conversion decay — > ir e + e~ using the complete data sample available at 
CMD-2. 
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2 Experiment 



The general purpose detector CMD-2 installed at the e + e collider VEPP-2M 



[ |15| has been described in detail elsewhere [ |i~6|j . 

It consists of a cylindrical drift chamber (DC) and double-layer multi- 
wire proportional Z-chamber, both also used for the trigger, and both inside 
a thin (0.38 Xq) superconducting solenoid with a field of IT. The momentum 

resolution of the DC is equal to a p /p = ^90 ■ (p(GeV)) 2 + 7% . The accuracy 
in the measurement of polar and azimuthal angles is oq = 1.5 • 10~ 2 and 
00 = 7 • 1CT 3 radians respectively. 

The barrel calorimeter with a thickness of 8.1Xq is placed outside the 
solenoid and consists of 892 Csl crystals. The energy resolution for photons 
is about 9% in the energy range from 50 to 600 MeV. The angular resolution 
is of the order of 0.02 radians. 

The end-cap calorimeter placed inside the solenoid consists of 680 BGO 
crystals. The thickness of the calorimeter for normally incident particles is 
equal to 13.4X . The energy and angular resolution varies from 8% to 4% 
and from 0.03 to 0.02 radians respectively for the photon energy in the range 
100 to 700 MeV. Both barrel and end-cap calorimeters cover a solid angle of 
0.92 x47r steradians. 

The experiment was performed in the <fi meson energy range (985-1060 
MeV). The integrated luminosity collected during the runs of 1993 (PHI93), 
1996 (PHI96) and 1998 (PHI98) was 1.5, 2.1 and 11.5 pb" 1 respectively, so 
that our analysis is based on the data sample corresponding to 15.1 pb _1 . 



3 Data analysis 
3.1 General approach 

A search for events of the decay <p — > 7r°e + e _ used the final state with two 
charged particles and two photons from the decay 7r° — > 77. For the process 
under study the invariant mass of the two photons M inv (77) should be equal 
to the 7r°-mass within the resolution. Event selection was performed using 
kinematical reconstruction taking into account energy-momentum conserva- 
tion. The number of detected events is given by: 

iV^Oe+e- = Nf B ((/>-> vr°e + e~) ■ B(ti 77) • £^o e + e - (1) 

where is the total number of the produced mesons and e^^o e + e - is the 
corresponding detection efficiency. 
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One of the most important backgrounds to our process comes from events 
of the decay ~ * 7r °7) 71-0 ~~ * 77 followed by the 7-quantum conversion at 
the wall of the beam pipe. Since the DC space resolution is not sufficient to 
separate events with the conversion at the beam pipe from those due to con- 
version decays at the interaction point, the contribution of this background 
was calculated from the simulation: 

iV^o 7 = N4 ■ 5(0 -> 7r° 7 ) ■ B(n° -> 77) • e^_ Mr o 7 , (2) 

where is the detection efficiency of the — > 7r°7 decay with the 7 

conversion in the material. 

Another significant source of background is the decay mode — > 7r + 7r~7r°. 
In this case the kinematical reconstruction assuming that charged tracks are 
electrons, gives M inv (77) about 150 MeV, i.e. only slightly different from the 
7r°-mass, and is not sufficient to separate this process from (|l|). Therefore, to 
suppress this background, the procedure of e/7r separation by the energy loss 
of charged particles in the calorimeter was developed (see the next subsection 
and [ |17| ) allowing determination of the expected number N^^+^-^o of such 
events. 

Other sources of background are the quantum electrodynamical (QED) 
process e + e~ — > e + e~77, and the decay mode — > 777 followed by the Dalitz 
decay 77^e + e~7 or the conversion of 7 from the 77 — > 77 decay. All these 
background processes can be separated from the decay under study by their 
wide and smooth distribution of M inv (77) . 

The total number of observed events is a sum of the three contributions 
above: 

N^Z^e+e- = A^->7r°e+e- + iV^_ >7r o 7 + N^^+^-^o. (3) 

The following selection criteria were used for the decay — > n°e + e~: 

• the impact parameter of the tracks d < 0.25 cm to reject events with 
the Kg L decay in DC; 

• the angle between two tracks Atp < 0.3 radians to suppress the process 
— > 7r + 7T~7r° (for conversion decays the angle between e + and e~ is 
small) ; 

• the number of photons N^ REC = 2; 

• the photon energy E™ ax < 460 MeV to suppress the process e + e _ — > 

e + e~77; 

• E™ in > 50 MeV for calorimeter noise suppression; 
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• \E™ ax — 362 1 > 15 MeV for suppression of the conversion decay 
7]—>-e + e~'-f. 

To determine the total number of the produced </> mesons, events of the 
process — > 777, 77 — ► 7r + 7r _ 7 have been used. Since this process also contains 
two charged particles and two photons in the final state, normalization to it 
allows to cancel some possible systematic errors. The number of events of 
this process is given by the following formula: 

N ri ^ 7T+n - J = ■ B((f) -> 777) • B{t] -> 7r + 7r~7) • e r ,_ +7r+7r - 7 . (4) 

For the normalization process the selection criteria are: 

• Atp < 2.5 to suppress K s — > tt + 7t~ events; 
. N* REC = 2; 

• Mi nv (77) > 250 MeV for suppression of <p — > 7r + 7r~7T° events; 

• d < 0.25 cm to reject events with the K SL decay in DC. 

Figure [l] shows the distribution of the invariant mass M inv (7r + 7r~7) for 
thus selected experimental events. A clear signal from the process — > 777, 
77 — > 7r + 7r _ 7 is observed at the 77-meson mass. It can be fit with a Gaussian 
whose variance is taken from simulation. The background mainly comes 
from events of the process (f> —> 777, 77 — > 7r + 7r~7r° with a lost photon (the 
wide signal in the 500 MeV region can be fit with the logarithmic Gaussian [ 
18]), the decay <fi — > 7r + 7r~7r° and the QED process e + e~ — > e + e~77 (another 



wide background fit with the Gaussian). 

Eliminating from the relations (1), (2) and (4), the following expres- 
sion for the branching ratio of the decay <fi — > ir°e + e~ can be obtained from 
(3): 

B((f) -> 7r°e + e-) = 

R(rh _> nVl ^5 ~^ 7T + 7r~7) g !Z - ! 7r+7 C 7 ^^ifie+e- ~ N <l>->n+n-x 

W V1) ' B(7r° -+ 77) ' e^rOe+e- ' A^ +7r - 7 (5) 



</>— >7r e+e — 

Here N^ n + n - n o is the expected number of the events of the decay 
7T + 7r~7r obtained from the procedure of e/n separation. 
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Figure 1: Invariant mass M inv (tt + 7t 7) for <fi — > 777, 77 — > tt + ti 7 
3.2 e / 7r-separat ion 

Since simulation can not completely reproduce interactions of charged par- 
ticles with the calorimeter, the procedure of e/7r-separation is based on the 
experimental data. Energy losses of and 7r =1= were studied using "clean" 
samples of events of the processes e + e~ — > e + e~7 and <fi — > 7r + 7r~7r°, 7r° — > 77 
respectively. The QED process e + e~ — *> e + e~7 with small angles between the 
tracks was chosen as the closest to <fi — > 7r°e + e~. 

From the distributions of E/P (the energy loss of a charged particle in 
the calorimeter divided by its momentum in DC) versus P shown in Fig. Q a) 
for pions from the — > 7r + 7r~7r° decay and in Fig. ||] b) for electrons from the 
e + e~ — > e + e~7 it can be seen that e/ir can be separated well. 

The whole region of momenta from to 500 MeV was divided to bins 
50 MeV wide and for each bin the E/P distributions were plotted. These 
spectra were fit by: for e - the logarithmic Gaussian function, for - 
a sum of two standard Gaussian functions. A probability to penetrate the 
detector without any energy loss and a probability to lose some energy in the 
calorimeter were defined as a ratio of the number of events with E — and 
with E > to the total number of events respectively. 

Then the probability density functions for each sort of the particles / e ± )7r ± 
depend on the momentum in DC and energy loss in the calorimeter. 

The probability for one track to be tt + or 7r~: 



(6) 



/tt+ + fe+ 



7T 



fn~ + fe 
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Figure 2: Distribution of E/P vs P: a) tt + for 

e + e~ — > e + e~7 



7r + 7r 7T°; b) e + for 



The probability for one track to be e + or e : 

W e+ = fe l , W e - = fe ~ (7) 

Jn+ + Je+ Jw- + Je- 

Similarly, one can construct the probability for both tracks in the event 
to be 7r + 7r~ or e + e~: 

. w w - + w e+ ■ w e - ' We+e ~ - • W v - + W e+ ■ W e - [ * } 

Finally, we can determine the detection efficiencies £ e J+ e - and £^.+ r 7r - which 
respectively give the probability to detect an e + e~ and 7r + 7r~ pair as an e + e~ 
one. They are calculated as the ratio of the number of events with W e + e - > 
Wo to the total number of sample events of the processes e + e~ — > e + e~7 and 
<f) — > 7r + 7r _ 7T° respectively. Here Wq is some boundary value chosen from the 
analysis of the W e + e - distribution. 



3.3 Selection of </> — > 7r e + e events 

The main selection criteria for the process under study were listed above. 

Figure § shows the scatter plot of the probability W e + e - for both tracks 
to be electrons versus the invariant mass of the photons Mi nv (77) for the 
PHI98 data after these cuts. Events with an e + e~-pair have W e + e - ~ 1 
while for those with pions W e + e - ~ 0. Events in the region W e + e - ~ and 
M inv (77) ~ 150 MeV come from the decay <fi — > 7r + 7r~7r°. To suppress it, the 
additional cut based on the e/7r-separation was applied: 
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Figure 3: Scatter plot of W e + e - versus M inv (77) 



• W e + e - > 0.5, the probability that an event has an e + e -pair. 

Figure |] shows the distribution of the invariant mass M inv (77) for the 
data after this selection. A clear signal is observed at the 7r°-meson mass 
which can be fit with a Gaussian (its variance is taken from simulation). 
The wide background can be also fit with a Gaussian. 

If instead the cut W e + e - > 0.5 one applies the complementary cut 

• W e + e - < 0.5, to enhance the probability that the event has a 7r + 7r~ 



the expected number N^^+^-^o can be determined. 

Detection efficiencies were determined from the Monte Carlo simulation 
(MC) [ [HJ] taking into account the efficiencies of e/ir separation £ e J+ e _,e'^+ 
and a correction for the reconstruction of events with a small open- 
ing angle, for brevity referred to as the small angle correction. The recon- 
struction of such events can not be completely reproduced by the simulation 
and was therefore specially studied using experimental events of the process 
<f) — > 7r + 7r~7T°, 7r°^e + e~7 which also have a small angle between an electron 
and positron. 

_ MC e/w _ MC e/7r 



pair, 



£ <j> — >7T+7r 7T° 



MC 

(f> — >7 
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Figure 4: The invariant mass M inv (77) for — > 7r°e + e , tt c 



77 



4 Results 

Figure |^ shows the energy dependence of the visible cross section of the 
process — > 7r°e + e~, 7r° — > 77 for the run PHI98. It is compatible with that 
expected for the meson, but a small non-resonant background can not be 
excluded. In fact, one should expect such a background coming from the 
"tail" of the uo — > 7r°e + e - decay. To determine the number of events coming 
from the meson decay, a fit of the visible cross section was performed which 
included a Breit-Wigner signal of the meson and a possible non-resonant 
background. The fit with the meson mass and width fixed at their world 
average values gave the following ratio of the non-resonant background cross 
section to that at the peak 

^ = 7.5 +15 7 l% 
04, 

compatible with the VMD estimate. From this value as well as from the 
integrated luminosity one can determine the corrected number of events for 
the process under study: 

N e / P 0p+p _ = 47.6 ±11.1. 

Similar analysis was performed for the data of PHI93 and PHI96 runs. 

Table [l] presents the final results of the data processing: the number of 
selected events, the expected number of events of — > 7i°e + e~ , — > 7r°7 and 
— > 7t + 7t~7t° determined by the described procedure, the detection efficien- 
cies and branching ratios of the — * n°e + e~ decay in various runs. Also 
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Figure 5: The visible cross section of — > 7r°e + e~, it — > 77 versus the total 
energy 



presented in the Table are the data for the normalization process. The dif- 
ferences of the efficiencies from run to run are caused by varying conditions 
of data taking. The total detected number of the — > n°e + e~ candidates is 
67.9 ± 12.3 with an expected background of 22.3 ± 1.4 events. The error in 
the estimated number of background events is determined by the procedure 
described above. 

As already noted, since the number of events of the process under study 
is normalized to that of the process <fi — > 777, 77 — > 7r + 7r~7 with a similar 
final state, some of the systematic uncertainties cancel (trigger, detection 
efficiencies etc.). 

The main sources of the remaining systematic errors are listed below, 
their magnitude given for the most statistically significant run PHI98: 

• A limited sample of simulated events used to determine detection effi- 
ciencies - 2.7% 

• Statistical errors of the parameters in the small angle correction - 3.9% 

• Parameters of the e/n separation procedure - 2.5% 

• The shape of the invariant mass distributions used in the fits - 11.9% 

• Inaccurate knowledge of the thickness of material in front of the DC 
(the beryllium beam pipe and the DC inner wall made of aluminized 
mylar) - 2.2% 
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Table V. Branching ratio of — > ir e + e decay 



Run 


PHI93 


PHI96 


PHI98 


Total 




6.3 ± 2.8 


14.1 ± 4.4 


47.6 ± 11.1 


67.9 ± 12.3 


±0.1 


±0.5 


±3.6 


±5.2 


C MC 0/ 


19.11 ±0.31 


26.74 ±0.37 


23.10 ±0.34 




£ 7rO e +e-> /0 


±2.18 


±2.51 


±1.82 






4.98 ±0.22 
±0.40 


7.30 ±0.27 
±0.59 


5.44 ±0.23 
±0.44 






99.5 ±20.1 
±7.1 


91.0 ± 10.4 
±5.6 


96.4 ± 2.8 
±4.7 




£„+„.-, /o 


92.4 ±8.3 


93.3 ±5.4 


89.9 ± 1.3 




9.1 ±2.9 


7.8 ± 1.0 


7.0 ±0.4 




Expected 


0.9 ±0.3 


2.2 ±0.3 


8.5 ±0.6 


11.7 ±0.7 


-^M^> — >7T°7 


±0.1 


±0.3 


±1.1 


±1.5 


Expected 


1.4 ±0.6 


1.7 ±0.5 


7.4 ± 1.0 


10.6 ± 1.3 




±0.1 


±0.1 


±0.1 


±0.1 


Total expected 
background 


2.3 ±0.7 


4.0 ±0.6 


16.0 ± 1.1 


22.3 ±1.4 


Expected 


3.9 ±2.9 


10.1 ±4.4 


31.6 ± 11.2 


45.6 ± 12.4 




12.34 ±0.25 


20.47 ±0.32 


20.75 ±0.32 




Nf] — >7V + TY~ 7 


126 ± 17 


362 ± 25 


1858 ± 58 


2346 ± 65 


± 6 


± 9 


±38 


±53 


B(0 -> 7r°e+e-), 


1.36 ± 1.13 


1.57 ±0.75 


1.10 ±0.40 


1.22 ±0.34 


10~ 5 


±0.22 


±0.23 


±0.19 


±0.21 



• The shape of the distributions used in the fits to determine the param- 
eters of the small angle correction - 6.8% 

• Errors of the world average values for the branching ratios of interme- 
diate decays from [ |9j - 4.9% 

• Dependence on the transition form factor model. The form factor in 
the generalized VMD (p ± p' mesons) was compared to that in the 
simple VMD (single p meson only). The contribution to the branching 
ratio is 9.5%. 

Three first sources of the error are of statistical nature, their values were 
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obtained for each run separately, so that they are uncorrelated. Their total 
contribution to the error is 5.4% and can be quadratically added to the 
statistical error. All others are obviously correlated from run to run and 
their total magnitude is 17.5%. 

From the last line of Table 1 it is clear that the values of the branching 
ratio obtained in different runs are consistent within the errors. To improve 
the accuracy, we can average them. The resulting value of the branching ratio 
for the — > 7r°e + e~ decay is (1.22 ± 0.34 ± 0.21) • 10~ 5 where the first error 
is statistical (including the uncorrelated systematic errors discussed above) 
and the second one is a systematic correlated error common for all runs. 



5 Conclusion 

For the first time the branching ratio of the conversion decay (p — > 7r°e + e~ 
has been determined by the CMD-2 detector at VEPP-2M: 

B(<p -> 7r°e + e-) = (1.22 ± 0.34 ± 0.21) • 10~ 5 . (9) 

This measurement is based on 68 selected candidates for the events of the 
process <fi — ► 7r°e + e~ with an expected background of 22 events. The obtained 
value agrees with the theoretical predictions (1.3 — 1.6) • 10~ 5 [ f|, ^(J, ^1| 
and with the experimental upper limit 1.2 • 10 -4 at 90% CL placed by the 
neutral detector ND at VEPP-2M using a data sample of 2.8 pb _1 [ 11 1. Our 
result supersedes the previous value of (0.85 ± 0.61 ± 0.12) ■ 10 -5 obtained 



by CMD-2 and based on a data sample corresponding to 1.5 pb only [ |22 



The applied procedure of event selection and particularly the cut on the 
angle between the tracks A\l/ < 0.3 selects events with a rather small q 2 so 
that it is practically impossible to study the momentum transfer dependence 
of the cross section. Therefore, one can conclude that the study of transi- 
tion form factors will become feasible when much larger data samples are 
collected. 

The authors are grateful to the staff of VEPP-2M for the excellent per- 
formance of the collider, to all engineers and technicians who participated in 
the design, commissioning and operation of CMD-2. We acknowledge useful 
comments of V.P. Druzhinin and stimulating discussions with R.A. Eichler. 
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